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Memory  Effects  In  Dynamical  Many-Body  Systems: 
The  Iaomneslc  (Constant-Memory)  Approximation 


A.  C.  Berl  and  Thomas  F.  George 
Department  of  Chemistry 
University  of  Rochester 
Rochester,  New  York  1^627 


Abstract 


Exact  numerical  solutions  of  a  generalized  master  equation  for  a  system 
involving  a  temporally  local izable  memory  kernel  and  a  constant  one  are 
presented.  Comparison  with  an  isomnesic  (constant-memory)  approximation,  with 
the  memory  kernel  reduced  to  the  sum  of  a  local  and  a  constant  term, gives 
excellent  agreement.  A  Markovian  version  of  the  problem  is  also  described  and 
found  to  be  inherently  ambiguous,  devoid  of  significant  physical  features  seen 
in  the  transient  regions,  but  leading  to  the  same  steady-state  as  the  isomnesic 
approximation. 


I.  Introduction 

The  dynamics  of  a  system  of  many  interacting  particles  is  often  treated  by 

partitioning  the  system  into  two  or  more  subsystems,  one  of  which, £5,  comprises 

the  degree(s)  of  freedom  of  primary  interest,  while  the  others,  js,  ,jt2  f  • • « i  3T € 

1 

reservoirs  to  which  ^  is  coupled.  No  observations  are  made  on  the  reservoirs, 
and  the  corresponding  degrees  of  freedom  are  projected  out  of  the  equations  of 
motion.  A  quantum  mechanical  generalized  master  equation  (GME)  results,  and 
invokation  of  the  Born  and  random-phase  approximations  generates  a  set  of 
coupled  Volterra-type  integrodiff erential  equations  (VTE) 


Ps(t) 


I 

S'*S 


fl 

I  df [k  ss(t-t')Ps,(t')  -  Kss,(t-t')  Ps(t’)] 


(1) 


where  P  (t)  is  a  diagonal  matrix  element  of  the  projected  density  operator  in 
the  subspacejo,  and  represents  the  probability  of  finding  Jo  in  an  eigenstate  |S> 
of  the  zero-order  sub-Hamiltonian  3V 

<3Cg|s>  -  6$|S>  .  (2) 

The  kernel  functions  Ksg, (t-f )  include  temporal  correlations  between  dynamical 
variables  of  the  reservoirs.  The  resulting  dependence  of  the  state  of  Jo  at  any 
time  t  upon  all  previous  times  t'  £  t,  the  memory  effect,  is  manifested  in  the 

3 

convolution  form  of  Kgs, (t-t1 ),  the  memory  kernels.  For  most  real  physical 
systems,  the  form  of  Kgs,  precludes  an  analytic  closed-form  solution,  and  often 

ij 

even  step-by-step  numerical  integration.  When  numerical  solutions  using 

5 

alternative  techniques  such  as  the  method  of  successive  approximations  are 


possible,  one  is  limited  to  small  times.  For  example,  for  kernels  with  some 
form  of  periodicity,  a  fairly  common  situation,  these  times  are  of  the  order  of 


a  few  periods,  which  seldom  correspond  to  macroscopic  measurement  times.  The 
solutions  are  therefore  of  little  use  in  suggesting  or  elucidating  physical 
mechanisms. 

The  need  for  numerical  solutions  has,  in  the  past,  been  circumvented  by 

6  7 

making  the  Markovian  approximation  (MA),  *  namely  that  the  system  is  amnesic 
(has  negligible  memory).  Thus  the  convolution-type  kernels  are  assumed  to  be 
extremeley  localized  temporally, 

KgS.U-t’)  -  Wss,  6 ( t — t * )  ,  (3) 

which  is  tantamount  to  assuming  that  the  probability  functions  are  very  slowly 

varying  during  the  (small)  period  of  time  in  which  Kss,(t)  changes 

substantially.  This  approximation  is  difficult  to  justify  in  general,  and  is  in 

fact  not  justified  if  a  locally  correct  solution  is  sought.  This  is  because 

Ps(t)  must  reflect  the  nature  of  Kss,(t)  itself,  and  if  the  latter  varies 

substantially  in  any  time  range,  so  must  Pg(t),  at  least  locally  in  time.  What 

makes  the  MA  useful  is  its  avowed  ability  to  provide  the  rough  behavior  of  Pg(t) 

over  long  time  periods.  Thus  on  a  coarse-grained  time  scale,  the  MA  is  expected 

to  provide  results  roughly  resembling  the  exact  results. 

Not  all  physical  systems  admit  of  the  representation  of  Kss,(t-t’)  shown  in 

Eq.  (3),  an  important  example  being  a  system  resonantly  pumped  by  a  continuous- 

7-9 

wave  laser  and  coupled  to  another  reservoir.  For  this  important  class  of 

physical  systems  the  kernels  can,  under  favorable  conditions  to  be  described 

later,  be  written  as  sums  of  terms  such  as  Eq.  (3)  and  additional  constants 

2 

independent  of  time.  The  resulting  equations  of  motion  are  very  different  from 
those  obtained  using  Eq.  (3).  The  differences  are  quantified  andhighlighted  in 
this  work,  and  it  is  shown  that  the  results  of  the  exactly-sol vable  equations 
with  kernels  including  the  constants  provide  an  excellent 


representation  of  the  solutions.  The  MA,  on  the  other  hand,  is  shown  to  lead  to 
probability  profiles  in  gross  disagreement  with  the  exact  results,  especially 
for  the  short  time  scale. 

The  format  of  this  paper  is  as  follows.  In  Section  II  the  basic  theory  is 
summarized,  and  the  results  of  calculations  are  presented  in  Section  III. 
Conclusions  are  given  in  Section  IV. 

II.  The  Isomnesic  (Constant-Memory)  Approximation 

A  system  of  great  current  interest  that  lends  itself  ideally  as  a  vehicle 

for  illustrating  the  various  approximations  for  the  memory  kernel  is  an  adatom 

adsorbed  on  a  solid  surface,  vibrationally  excited  by  an  infrared  continuous- 

2  7-9 

wave  laser,  and  coupled  to  phonon  modes  of  the  solid.  ’  The  system  9Sis  the 

adbond,  viz,  the  bond  resulting  from  the  interaction  of  the  adatom  with  the 

ohonon-averaged  motion  of  the  lattice  atoms;10  the  phonons  and  the  laser  are  two 

independent,  noninteracting  reservoirs^  andJJJ^.  Quantum  theoretical 

treatments  of  this  class  of  systems  lead  quite  generally  to  a  GME  and  a 

corresponding  VTE,  Eq.  (1),  with  P<,(t)  representing  the  probability  of  finding 

the  adbond  in  the  eigenstate  |S>  of  the  zero-order  adbond  Hamiltonian  3^.  The 

2 

kernels  KSg,(t-t')  incorporate  the  influence  of  the  phonons  and  the  laser.  We 

have  partitioned  the  total  Hamiltonian  of  the  system  in  such  a  way  that  the 

perturbations  which  generate  the  kernels  are  traceless  in  their  individual 

2  1 0 

reservoir  subspaces,  ’  As  a  direct  consequence,  the  total  kernels  separate 
into  individual  subspace  kernels, 

KSS'^  "  KSS'^  +  KSS'  (A) 

where  Kg^(t)  is  due  to  the  phonon  field  alone  and  Kgg)(t)  is  due  to  the  laser 


radiation  field  alone.  Using  a  one-dimensional  model  for  the  solid, 


we  have  obtained  the  kernels  in  closed  form  and  solved  them  numerically  by  an 

5 

iterative  technique  for  times  of  the  order  of  a  few  "periods"  of  Kss,(t). 

These  periods  are  dictated  by  the  characteristic  frequencies  of  the  phonons  and 
the  adbond  (i.e.,  the  Debye  frequency  ui-  and  the  transition  frequencies  w  for 

u  oo 

bound  states  of  the  HamiltonianSCg) ,  and  the  laser  frequency  all  of  which  are 
12  1 1) 

-10  -  10  Hz  (1-100  THz).  The  periods  are  therefore  very  small,  being 

-12  -I1) 

10  -  10  s  (1-.01  ps);  consequently,  the  solutions  of  the  VTE  can  be 

obtained  numerically  for  times  no  larger  than  about  10  ps,  which  is  certainly 

insufficient  for  the  description  of  a  physical  or  chemical  process.  In  order  to 

obtain  solutions  for  longer  times,  an  analytic  solution  must  be  sought. 

We  begin  with  the  assertion  that  the  physically  important  situation  is  one 

where  the  laser  frequency  and  the  Debye  frequency  equal  different  transition 

2  (r ) 

frequencies  wss,  of  the  adbond.  For  this  case  Kss,(t)  becomes  a  constant  , 

(from  the  term  with  |^ss,|  *  P*us  oscillating  terms  which  we  ignore  in  the 

sense  of  the  rotating-wave  approximation.^  The  phonon  term,  K^)(t),  has  a 
complicated  behavior  typified  by  K^^(t)  in  Fig.  1,  but  we  note  that  it  has 
appreciable  magnitude  for  times  -  2  ps,  and  almost  vanishes  subsequently.  If 
the  mesh  size  of  our  theoretical  experiment  is  no  less  than  2  ps,  we  can  treat 
K^^(t)  as  a  delta  function,  so  that,  from  Eq.  (H), 

KSS,U>  -  0SS,  5<t>  •  kss,  .  (5) 

where 

f 

°SS'  ■  \  at  Kss’  (t)  •  (6) 

'o 

In  Eq.  (5),  the  first  term,  originating  from  the  coupling  of  the  adbond  to  the 
phonons,  is  the  amnesic  or  Markovian  term,  and  would  lead,  in  the  absence  of  the 
laser  term,  to  the  elimination  of  all  memory  effects.  The  second  term,  due  to 


o 


the  resonant  laser  coupling,  is  the  isomnesic  non-Markovian  term,  and  represents 
a  constant  memory  effect  all  the  way  back  to  t-0. 

The  subsequent  treatment  is  simplified  by  defining  "diagonal"  elements  Uss 
and  ksg,  viz. 


-  I  Q. 

S'*S 


SS1 


(7) 


k 


SS 


and  the  notation 


-I 

S'*S 


kS’S 


(8) 


n  e  lQss,l  ,  (9a) 

k  *  ikss,}  ,  (9b) 

K  (t)  s  {Kss,(t)}  ,  (9c) 

W  £  iwss,}  ,  (9d) 

-  *  ^6SS’^  ’ 

P  (t)  s  {Ps(t)l  .  (9f) 

The  exact  VTE  may  be  written  as 


ft 

P  (t)  «  I  dt'  K  (t-t’ )*£(t')  ;  (10) 

o 

in  the  isomnesic  approximation  (IA)  it  becomes 

P  (t)  -  k  •  I  dt'  P  (t')  +  O-P(t)  .  (11) 

^0 


Formally  this  can  be  solved  without  any  further  approximations.  Differentiating 
both  sides  of  Eq.  (11),  we  get  the  second-order  differential  equation 


P(t)  -  k’P(t)  +  0«P(t)  , 


(12) 
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whose  Laplace  transform, 

M(s)-f(s)  ■  T(s)  ,  (13) 

where 

M(s)  s  k  ♦  SO  -  s2l  ,  ( 1 1») 

T(s)  E  -  sP(0)  ,  (15) 

and£(s)  is  the  Laplace  transform  of  P(t),  is  easily  inverted  using  the 
Heaviside  expansion  theorem.^1 

While  this  effectively  eliminates  the  non-Markovian  bottleneck  without 

recourse  to  Eq.  (3),  it  is  instructive  to  generate  solutions  of  a  Markovian 

version  of  this  class  of  problems.  To  do  this,  we  have  to  make  provisions  for 
( r ) 

localizing  K^,  in  time.  This  brings  to  the  fore  the  pervasive  issue  of  time 

scales  in  dynamical  many-body  systems,  our  specific  question  being  the  relative 

degree  of  localization  ascribable  to  30(3  KSS’’  the  laser  radiation  is 

localized  as  a  pulse,  it  can  still  extend  over  a  period  of  time  much  longer  than 

the  extent  of  K^g)(t).  In  other  words,  the  delta  function  representing  Kgj^(t) 

( r ) 

may  be  quite  different  from  the  delta  function  representing  Kgg,(t),  and  the 
former  may  in  fact  not  really  be  considered  as  localized  with  respect  to  the 
latter.  With  these  provisos  in  mind,  we  can  proceed  to  write  the  Markovian 
version  of  Eqs.  (5),  (10),  (12)  and  03): 


K(t)  -  (n  ♦  f)  «(t)  , 

(16) 

P(t)  -  w  •  P  (t)  , 

07) 

J(s)  •  fi(s)  •  u  , 

(’8) 

where 


f  -  k  t  , 


(19) 
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w  -  a  *  f  ,  (20) 

m  m  m 

J(s)  -  (Q  ♦  f  -  sA)  ,  (21) 

m  *  «  «  * 

U  -  -  P(0)  ,  (22) 

and  t  represents  the  actual  extent  of  the  laser  signal  or  a  variable  time 
parameter.  Equation  (17)  is  the  well-known  Pauli  master  equation  (PME).1’^ 

To  summarize  the  formulation,  Eq.  (10)  is  the  exact  GME  with  a  convolution 
type  kernel  which  makes  solutions  dificult  for  experimental  times.  If  the 
kernel  can  be  approximated  as  the  sum  of  a  very  sharply  peaked  function  due  to 
the  phonons  plus  a  constant  due  to  the  laser  (the  isomnesic  approximation),  we 
get  the  simplified  GME,  Eq.  (11),  which  is  still  not  Markovian  but  has  a 
"constant  memory"  (isomnesic-)  term.  The  resulting  second-order  differential 
equation,  Eq.  (12),  can  be  solved  analytically  for  arbitrarily  long  times.  If, 
on  the  other  hand,  memory  effects  are  assumed  to  be  negligible,  i.e.,  with  a 
kernel  highly  localized  in  time,  Eq.  (3).  we  get  the  Markovian  limit,  namely  the 
PME,  Eq.  (17),  which  is  a  first-order  differential  equation  and  can  also  be 
solved  analytically  for  arbitrarily  long  times.  The  nature  of  the  solutions 
depends,  of  course,  on  the  time  paramter  t  in  Eq.  (19). 

III.  Results  of  the  Exact,  Isomnesic  and  Markovian  Approximations 

For  a  linear-chain  model  of  the  solid,  the  exact  kernel  function  due  to 
phonons  alone,  K^,  involves  an  expansion  over  the  number  of  phonons  n  and  the 
number  of  lattice  atoms  l.  For  the  specific  system  chosen,  the  effective 
surface  potential  generates  five  bound  states;  we  have  chosen  |  -  ujd  and 
|«  |  »  <*>L  [u>gs »  '  (&s-£s,)/tf;  see  Eq.  (2)].  The  phonon  kernels,  e.g.,  K^(t) 

in  Fig.  1,  fall  off  quite  rapidly  in  time  and  arise  almost  entirely  from  the 


(n-1  ,  1-1 )  term. 
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At  issue  is  the  validity  of  the  isomnesic  and  Markovian  approximations.  We 
can  examine  two  aspects  of  the  results  ,  viz.  the  details  of  the  early  transient 
behavior,  and  the  long-term  steady  state.  For  the  Markovian  case,  we  can 
further  examine  the  effects  of  varying  t  in  Eq.  (19).  The  transient  behavior 
can  be  compared  with  exact  numerical  solutions  of  Eq.  (10),  but  the  steady-state 
results  of  the  MA  and  the  IA  can  only  be  judged  by  extrapolation. 

While  the  probability  profiles  P(t)  provide  sufficient  comparison,  the 
physical  contents  of  our  results  are  better  displayed  by  the  average  adbond 
energy  given  by 

£(t)  =  l  Pc(t)6Q  ,  (23) 

S  b  b 

where  6S  is  the  energy  eigenvalue  corresponding  to  the  level  |S>  of  the 
effective  adbond  potential.  This  provides  a  simple  measure  of  the  chemical 
direction  of  the  adbond  and  a  well-defined  connection  between  this  exact  theory 
and  phenomenological  models  of  such  processes. 

The  results  obtained  by  solving  Eqs.  (10),  (12)  and  (17)  for  a  variety  of 
situations  of  physical  interest  are  very  consistent,  and  we  present  a  typical 
set  of  curves  in  Fig.  2  for  early  times,  t  S  5  ps.  The  exact  results  are  rich 
in  detail  and,  to  a  large  degree,  follow  the  variations  in  K(t).  We  shall  not 
dwell  on  the  physical  contents  except  to  mention  the  conspicuous  absence  of  the 
simple  exponential  approach  to  equilibrium  postulated  in  relaxation  time 
theories.  The  results  of  the  IA,  while  exhibiting  the  expected  loss  of  the  fine 
detail,  are  seen  to  chart  out  the  average  behavior  of  £(t)  with  great  accuracy. 
Thus,  for  example,  the  crossing  point  of  P1 (t)  and  P^(t)  is  reproduced  at  t  ■ 

1.7  ps  as  compared  to  the  exact  crossing  point  at  2.2  ps,  and  the  difference 
P^t)  -  P^(t)  is  seen  to  follow  the  exact  difference  with  a  similar  slight 


10 


temporal  retardation.  The  gross  oscillations  of  P(t)  are  also  reproduced 

faithfully.  Consequently,  we  feel  confident  in  our  prediction  that  the  steady- 

state  results  of  the  IA  will  conform  to  the  exact  results  very  accurately.  We 

wish  to  emphasize  the  significance  of  this.  The  IA  makes  possible  the  com pi ete 

solution  of  a  GME  for  arbitrarily  long  times,  when  the  solution  of  the  exact  GME 

is  essentially  impossible.  The  solution  of  the  GME  in  the  IA  is  no  more 

difficult  than  that  of  the  PME  obtained  by  invoking  the  MA  which,  however,  gives 

results  which  differ  radically  from  the  exact  results  in  the  transient  region 

for  a  wide  range  of  values  of  t.  Figure  3  illustrates  the  variations  of  P  ( t ) 

-5  5 

with  t  in  the  MA  (10  ps  i  x  i  10  ps).  The  probability  profiles  are  seen  to 

be  insensitive  to  increases  in  i  beyond  -  10  ps,  which  corresponds  to  almost 

-2 

instantaneous  saturation.  The  apparent  insensitivity  of  P(t)  for  t  <  10  ps 

seen  in  Fig.  3  is,  of  course,  misleading,  since  in  this  region  the  P<,(t)  vary 

with  time  too  slowly  for  the  change  to  be  evident  in  the  20  ps  time  range  of  the 

plot.  Drawn  for  2  ns  in  Fig.  4  [we  present  only  P^(t)],  the  differences  become 

-2 

more  pronounced  for  t  <  10  ps,  and  the  strong  dependence  of  the  apparent 
"equilibration  time"  on  t  emerges  clearly.  We  see  no  unambiguous  way  of 
choosing  an  appropriate  value  for  t.  For  systems  which  are  known  or  expected  to 
equilibrate  very  quickly,  one  could  make  an  ad  hoc  semiempirical  choice  of  t  to 
solve  a  PME  for  the  steady  state.  However,  if  one  is  interested  in  the  details 
of  the  short-time  behavior  of  the  system,  one  must  solve  the  GME  using  the  IA. 

For  comparison  of  the  MA  with  the  IA  we  have  chosen  t  -  1  ps.  The  resulting 
profiles  for  the  MA  are  shown  in  Fig.  5,  along  with  our  IA  results  for  20  ps. 

The  complete  absence  of  oscillatory  behavior  representing  beats  of  the  system 
and  of  population  crossings,  and  the  presence  of  predominantly  exponential 
variation  of  P(t)  representative  of  a  "relaxation  time"  exhibited  by  the  MA 
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results  of  Fig.  5  are  not  necessarily  general  features  of  the  solutions  of  a 

PME.  Since  the  matrices  M  and  J  are,  in  general,  asymmetric,  both  the  GME  and 

the  PME  can  generate  oscillatory  solutions  P(t),  contradicting  earlier 
1 2 

contentions.  However,  the  total  lack  of  resemblance  between  the  early  PME 
results  and  the  exact  ones  for  the  present  case  is  alarming  and  raises  urgent 
and  serious  questions  about  the  validity  of  the  Markovian  approximation,  even 
for  a  rough  qualitatiave  estimate  of  the  early  transient  behavior  of  many-body 
systems  with  a  nonrandom  force.  At  best,  these  results  suggest  the  need  for 
great  caution  in  using  the  MA.  In  retrospect  this  is  not  surprising,  because 
the  representation  of  a  time- independent  cw  laser  beam  by  an  ultrashort  pulse  is 
tenuous,  whereas  that  of  the  very-short-range  (in  time)  phonon  correlation  by  a 
delta  function  is  quite  reasonable. 

The  results  for£(t),  shown  in  Fig.  6  for  the  exact  GME  and  the  IA  and  in 
Fig.  7  for  the  MA  and  the  IA  further  support  these  contentions.  The  IA 
generates  an  energy  profile  which,  except  for  a  slight  temporal  retardation  f  is 
a  low-resolution  representation  of  the  exact  result.  The  significant  increase 
in  the  bond  energy  at  the  expense  of  both  the  laser  and  the  solid  is  clearly 
evident  in  the  IA  results  of  Figs.  2  and  6.  The  results  of  the  MA  are 
qualitatively  different  from  the  exact  and  IA  results,  and  describe  a  monotonic 
approach  to  the  steady  state.  The  exact  and  IA  results  show  actual  periodic 
energy  transfer  out  of  the  adbond  followed  by  successive  build  up  due  to  overall 
energy  Inflow.  The  apparent  correspondence  seen  in  Fig.  7  between  the  magnitude 
of  6(t)  in  the  MA  and  the  "base  line"  about  which  the  results  of  the  exact  GME 
and  the  IA  oscillate  is  entirely  due  to  our  choice  of  t  ■  1  ps.  That  no  such 
correspondence  can  be  expected  in  general  is  evident  from  Fig.  8  showing  the 
variations  of£(t)  with  t.  The  steady-state  (t+»)  values  of  P(t)  and  £(t)  in 


1 2 


the  MA  have  been  found  to  be  the  same  as  those  obtained  with  the  IA  for  all 
cases  considered  so  far  (when  the  values  of  t  chosen  lead  to  steady  states  in 
the  time  ranges  considered). 

IV.  Conclusion 

Our  results  for  the  time  evolution  of  a  dynamical  many-body  system  subject 
to  a  variety  of  forces  point  out  the  inadequacies  of  the  Markovian  approximation 
when  the  resulting  memory  kernel  is  nonlocal  in  time.  As  an  alternative,  we 
present  the  isomnesic  approximation  which  is  characterized  by  a  constant  memory; 
this  eliminates  the  non-Markovian  bottleneck  associated  with  the  convolution- 
type  memory  kernel  in  the  Volterra-type  integrodif ferential  equation  resulting 
from  the  generalized  master  equation.  Closed  solutions  of  the  latter  are  in 
excellent  agreement  with  the  exact  results,  in  contrast  to  those  of  the 
Markovian  approximation. 
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Fig.  1.  The  phonon  part  KQ1  (t)  of  the  memory  kernel. 


Fig.  2.  Probability  profiles  Pg(t)  for  levels  S-  0,  1  and  3  of  the  adbond. 

Solid  lines:  Exact  numerical  solutions  of  the  GME;  broken  lines:  closed- 
form  solutions  of  the  GME  with  the  IA. 

Fig.  3.  P„(t)  for  the  system  of  Fig.  2  with  the  MA  for  a  ten-decade  range  of  t. 
(S)-P0(t);  (B)-P1(t);  <C)-P3(t). 

Fig.  4.  P^(t)  with  the  MA  for  2  ns. 

Fig.  5.  Comparison  of  Pg(t)  with  the  IA  (solid  lines)  and  the  MA  ( x«1  ps;  broken 
lines)  for  S«0,  1  and  3 

Fig.  6.  Average  adbond  energy  £(t)  obtained  via  the  numerical  solution  of  the 
exact  GME  (solid  line)  and  closed-form  solution  with  the  IA  (broken 
line). 

Fig.  7.  6(0  with  the  IA  (solid  line)  and  the  MA  (t«1  ps;  broken  line). 

Fig.  8.  Variation  of  6(0  with  t  in  the  MA. 
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